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We approximate the solution of some linear systems of SDEs driven by a fractional Brownian 
motion B H with Hurst parameter H £ (|, 1) in the Wick-Ito sense, including a geometric frac- 
tional Brownian motion. To this end, we apply a Donsker-type approximation of the fractional 
Brownian motion by disturbed binary random walks due to Sottinen. Moreover, we replace the 
rather complicated Wick products by their discrete counterpart, acting on the binary variables, 
in the corresponding systems of Wick difference equations. As the solutions of the SDEs admit 
series representations in terms of Wick powers, a key to the proof of our Euler scheme is an 
approximation of the Hermite recursion formula for the Wick powers of B H . 

Keywords: discrete Wick calculus; fractional Brownian motion; weak convergence; Wick-Ito 
integral 

1. Introduction 

A fractional Brownian motion B H with Hurst parameter H € (0, 1) is a continuous zero- 
mean Gaussian process in K with stationary increments and covariance function 

E[B«B?] = l(\t\ 2H + \sf H -\t- S \ 2H ). 

The process B 1 ^ 2 is a standard Brownian motion, but a fractional Brownian motion is not 
a semimartingale for H 7^ |. In this paper, we restrict ourselves to the case H > 1/2, in 
which the corresponding fractional Gaussian noise (B^ +1 — B^) n& -^ exhibits long-range 
dependence. 

In recent years, a lively interest in integration theory with respect to fractional Brow- 
nian motion has emerged (see, e.g., the monographs by Mishura or Biagini et al. [4, 15]). 
One of the extensions of the ltd integral beyond semimartingales is the fractional Wick- 
Ito integral. It is based on the Wick product o, which has its origin as a renormalization 
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operator in quantum physics. In probability theory, the Wick product with ordinary dif- 
ferentiation rule imitates the situation of ordinary multiplication with Ito differentiation 
rule (cf. Holdcn et al. [11]). Actually, this makes it a natural tool to apply for extending 
the Ito integral. 

We first consider the fractional Doleans-Dade SDE dS t = S t d° B^ 1 , Sq = 1, in terms 
of the fractional Wick-Ito integral. The well-known solution, exp(i?f r — ^t 2H ), is the 
geometric fractional Brownian motion, also known as the Wick exponential of fractional 
Brownian motion. Note that the Wick exponential has expectation equal to one and can 
therefore be interpreted as a multiplicative noise. Moreover, the ordinary exponential can 
be obtained from the Wick exponential by a deterministic scaling. Neither process is a 
semimartingale for H ^ h. The name "Wick exponential" is justified by the fact that it 
exhibits a power series expansion with Wick powers (i?f r ) ofe instead of ordinary powers. 

More generally, we consider a linear system of SDEs, 

dX t = {A x X t + A 2 Y t ) d*Sf , X = x , 

(1) 

dY t = {BxX t + B 2 Y t ) d*B?, Y = y . 

One can obtain Wick power series expansions for the solution of this system, too. Our goal 
is to approximate these Wick analytic functionals of a fractional Brownian motion. To this 
end, we require an approximation of a fractional Brownian motion and an approximation 
of the Wick product. 

There are several ways to approximate a fractional Brownian motion. One of the first 
approximations was given by Taqqu [22] in terms of stationary Gaussian sequences. We 
refer to Mishura [15], Section 1.15.3, for further approaches to weak convergence to a 
fractional Brownian motion. Sottincn constructed a simple approximation of a fractional 
Brownian motion on an interval for H > ^ by sums of square- integrable random variables 
in [21]. He used the Wiener integral representation of a fractional Brownian motion 
on an interval, B^ 1 = J Q zii{t,s)dB s , for a suitable deterministic kernel zn(t,s), due 
to Molchan and Golosov, and Norros et al. [16-18]. For this purpose, he combined a 
pointwise approximation of the kernel Zii(t,s) with Donsker's theorem. This approach 
was extended by Nieminen [19] to weak convergence of perturbed martingale differences 
to fractional Brownian motion. We shall utilize Sottinen's approximation with binary 
random variables throughout this paper. 

The main problem of applying the Wick product on random variables with continuous 
distributions is that it is not a pointwise operation. Thus, an explicit computation of 
the Wick-Ito integral is only possible in rare special cases. But this is precisely the 
advantage of the binary random walks. In such a purely discrete setup, we apply the 
discrete counterpart of the Wick product as introduced in Holden et al. [10]. Starting 
from the binary random walk, one can build up a discrete Wiener space, and the discrete 
Wick product depends on this discretization. This Wiener chaos gives the analogy to the 
continuous Wick products. For a survey on discrete Wiener chaos, we refer to Gzyl [9]. 
However, we will introduce the discrete Wick product in a self-contained way in Section 
3. 
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We can now formulate a weak Euler scheme of the linear system of SDEs (1) in the 
Wick-Ito sense, 

X? = XP_, + (A 1 XP_ 1 + A 2 Y l n _ 1 ) o n (Bffi - 5ji n 1)/n ), 
X£ = x , l = l,...,n, 

(2) 

Y? = Y^ + {B l Xf_ 1 + B 2 Y l " 1 )o n (Bffi B*? 1)/n ), 
Y o=Vo, l = l,...,n, 

where o„ is the discrete Wick product and (By'™ — B^^.,) are the increments of the 
disturbed binary random walk. As a main result, we show that the piecewise constant in- 
terpolation of the solution of (2) converges weakly in the Skorokhod space to the solution 
of (1). This is the first rigorous convergence result connecting discrete and continuous 
Wick calculus of which we are aware. As a special case, (2) contains the Wick difference 
equation 

XT = X?_ 1 +X?_ 1 <> n (B? / '»-B%? 1)/n ), X% = 1, 1 = 1,..., n. (3) 

As a consequence, the piecewise constant interpolation of (3) converges weakly to a 
geometric fractional Brownian motion, the solution of the fractional Dolcans-Dadc SDE. 
This was conjectured by Bender and Elliott [3] in their study of the Wick fractional 
Black-Scholes market. 

In [21], Sottincn considered the corresponding difference equation in the pathwise sense, 
that is, with ordinary multiplication instead of the discrete Wick product: 

X? = Xf-i + X?-i(B?fi - £(?.!")/„)' *o = 1> I = 1, • • ■ , »• (4) 
The solution is explicitly given by the multiplicative expression 

*P = ri(l + -<-\)/J)- (5) 

j'=i 

By the logarithmic transform of ordinary products into sums and a Taylor expansion, 
one obtains an additive expression for ln(A"™) which converges weakly to a fractional 
Brownian motion. In this way, Sottinen proved the convergence of X to the ordinary 
exponential of a fractional Brownian motion [21], Theorem 3. This approach fails for the 
solution of (3) since, in a product representation, analogous to (5), the discrete Wick 
product o n appears instead of ordinary multiplication. There is, however, no straightfor- 
ward way to transform discrete Wick products into sums by application of a continuous 
functional. 

However, the solution of (2) exhibits an expression which is closely related to a discrete 
Wick power series representation. Therefore, the convergence can be initiated explicitly 
for the Wick powers of a fractional Brownian motion, which fulfill the Hcrmite recursion 
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formula. We obtain a discrete analog to this recursion formula for discrete Wick powers 
of disturbed binary random walks. Actually, the weak convergence of these discrete Wick 
powers is the key to the proof for our Euler scheme. 

The paper is organized as follows. In Section 2, we give some preliminaries on the 
Wick-Ito integral with respect to a fractional Brownian motion and introduce the Wick 
exponential and other Wick analytic functionals. We then define the approximating se- 
quences and state the main results in Section 3. Section 4 is devoted to some L 2 - estimates 
of the approximating sequences. We prove convergence in finite-dimensional distributions 
in Section 5 and tightness in Section 6. 

2. Wick exponential and Wick analytic functionals 

In this section, we introduce the Wick product and the Wick-Ito integral, and describe 
the Hermite recursion formula for Wick powers of a zero-mean Gaussian random variable. 
We then obtain the Wick power series expansions for the solutions of SDEs (1). 

We consider a geometric fractional Brownian motion or the so-called Wick exponential 
of a fractional Brownian motion exp(B^ — ^t 2H ). For H = 5, this is exactly a geometric 
Brownian motion, also known as the stochastic exponential of a standard Brownian 
motion. For all H <G (0,1) and t > 0, it holds that t 2H = E[(Sf ) 2 ] and thus the Wick 
exponential generalizes the stochastic exponential. It is well known that exp(£> t — ^i) 
solves the Doleans-Dade equation 

dS t = S t dB u 5o = l, 

where the integral is an ordinary ltd integral. Actually, the Wick exponential of fractional 
Brownian motion solves the corresponding fractional Doleans-Dade equation 

dS t = S t d°B?, Sq = 1, 

in terms of a fractional Wick-Ito integral (cf. Mishura [15], Theorem 3.3.2). We want to 
approximate solutions of similar SDEs. 

Let $ and \1/ be two zero-mean Gaussian random variables. The Wick exponential is 
then defined as 

cxp <> ($) :=cxp($- iE[|$| 2 ]). 
For a standard Brownian motion (B t )t>o an d s < t < u, it holds that 

cxp°(B u - B t ) exp°(B t - B„) = cxp°(B u - B s ). 

Forcing this renormalization property to hold for all, possibly correlated, $ and VP, leads 
to the definition of the Wick product o of two Wick exponentials: 



exp°($) oexp°(^) :=exp°($ + *). 
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The Wick product can be extended to larger classes of random variables by density 
arguments (cf. [2, 6, 20]). For a general introduction to the Wick product, we refer to the 
monographs by Kuo and Holden et al. [11, 14] and Hu and Yan [13]. Note that the Wick 
product is not a pointwise operation. If we suppose that $ ~7V(0, <r), then we have, by 
definition, $°° = 1 , $ o1 = <f> and the recursion 



on+l 



Observe that it holds that 



— exp°(x4)) 
ax 



x=0 



— exp(x-$--E[|a;$| 2 ] 



2=0 



($ - xa ) exp ( x$ - -E[|x$| 



Suppose we have 



dw k 



exp°(w$) 



u;=0 
d2 _ dz 



for all positive integers k <n. Then, with z = w + x, ^ = 4^ = 1, we get 



<j)0("+l) 



dw 

d™ d 

du; n dec 



exp°(u;$) 



tu=0 



o — exp°(a;$) 
ax 



exp°((w + x)$) 



UJ=0,2=0 



a:=0 
d n+l 

dz^ 



r exp°(2:$) 



2 = 



We now obtain, by differentiation and the Leibniz rule, the following Wick recursion 
formula: 



dui r 



(<f> - wcr^) exp ( w$ - iE[|w<f>| 2 ] 



w=0 



d" 

($ - w 2 )- exp( iu$ - -E[|x$| 2 ] 

dw™ 

+ n(-tr cxp( w$- ^E[|w$| 2 ] 



iu=0 



(6) 



Define the Hermite polynomial of degree n £ N with parameter a 1 as 
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The series expansion 

ex J x -laA=jr ) ±h n Ax) (7) 



n=0 



then holds true. The first Hcrmite polynomials arc h a 2 (x) — 1, h^ 2 (x) = x. By the Leibniz 
rule, we obtain the Hermite recursion formula 

h n +\x) = xh^(x) - na 2 h n a2 -\x). (8) 

By the equivalent first terms and recursions (6) and (8), we can conclude that for any 
Gaussian random variable $ ~A/"(0, cr) and all n € N, we have 

= (9) 

By (7), we additionally have 

exp*($) = V -r**". (10) 

The fractional Wick-Ito integral, introduced by Duncan e£ aZ. [6] , is an extension of the 
Ito integral beyond the semimartingale framework. There are several approaches to the 
fractional Wick-Ito integral. Essentially, these approaches are via white noise theory, as in 
Elliott and von dcr Hock [7], and Hu and 0kscndal [12], by Malliavin calculus in Alos et al. 
[1], or by an S-transform approach in Bender [2]. In contrast to the forward integral, the 
fractional Wick-Ito integral has zero mean in general. This is the crucial property for an 
additive noise. The Wick-Ito integral is based on the Wick product. For a sufficiently good 
process (X s ) sG [ ^, the fractional Wick-Ito integral with respect to fractional Brownian 
motion (B^ )m,t] can be easily defined by Wick-Riemann sums (cf. Duncan et al. [6] or 
Mishura [15], Theorem 2.3.10). If we suppose that 7r„ = {0 = to < t\ < ■ ■ ■ < t n = t} with 
max^g^ |tj — ti—i | — > for n — > oo, then 

f X s d»Sf := lim J2 X U-> o (B* - B^), 

if the Wick products and the i 2 (51)-limit exist. For more information on Wick-Ito in- 
tegrals with respect to fractional Brownian motion, we refer to Mishura [15], Chapter 
2. 

By the fractional Ito formula (cf. [2], Theorem 5.3 or [4], Theorem 3.7.2), we have 

d(B*y k = k(B*)o h - 1 d*B* 1 , (B^r k = l {k=0} . (11) 
For the Wick exponential 

exp*(Bf) = ^-(Sfr fc , (12) 

fc=0 
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we obtain, by summing up the identity (11), the fractional Doleans-Dade equation, 

dS t = S t d°B?, S = l. (13) 
For any analytic function F(x) = X^fcLo IT 2 ^' we define the Wick version as 

oo 

F°($)=^^$ ofc . 

From the recursive system of SDEs (11), we obtain SDEs for other Wick analytic Junc- 
tionals of a fractional Brownian motion 

oo 

Recall the linear system of SDEs (1), 

dX t = (A 1 X t + A 2 Y t )d*Bf I , X = x , 
dY t = {B 1 X t + B 2 Y t ) d»Bf, Y =y a . 

The coefficients of the solution, 

**=E^) ofc ' ^=E^) ofc ' ( 14 ) 

k=0 ' k=0 

can be obtained recursively via (11) to be 

a =x , b Q = y , a k = Axdk-i + A 2 b k -i, b k = Bia fc _i + B 2 b k ^ 1 . 

Note that it holds that \a k \, \b k \ < C k for a C G R+. This is according to the recursive 
derivation of the coefficients and it ensures that the Wick analytic functionals X t and Y t 
are square- integrable (cf. the proof of Proposition 6). 



3. The approximation results 

Here, we present the approximating sequences and discuss the main results. More pre- 
cisely, we introduce the Donsker-type approximation of a fractional Brownian motion 
and the discrete Wick product, and obtain Wick difference equations, which correspond 
to the SDEs. We shall work with a fractional Brownian motion on the interval [0, 1], but 
all results extend to any compact interval [0,T]. 

Wc first consider the following kernel representation of a fractional Brownian motion 
on the interval [0, 1], based on works by Molchan and Golosov [16, 17], 

Bf= f z H (t,s)dB s . (15) 
Jo 
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For H > | , the deterministic kernel takes the form 



Z H (t, s) = l {t > s} c H (h-£) s^- h J* u H -^- (u s) H -^ dt 



(16) 



with the constant 



ch 



I 2HT(3/2-H) 
T{H + l/2)r(2-2ii") 



where T is the Gamma function (Norros et al. [18] or Nualart [20], Section 5.1.3). In 
order to simplify the notation, we think of Ha (^,1) as fixed from now on and omit 
the subscript H in the notation of the kernel. For an introduction to some elementary 
properties of fractional Brownian motion, we refer to Nualart [20], Chapter 5, Mishura 
[15] or Biagini et al. [4]. 

We apply Sottincn's approximation of a fractional Brownian motion by disturbed binary 
random walks. Suppose (f2, J 7 , P) is a probability space and, for all n £ N and i = 1, . . . ,n, 
we have independent and identically distributed symmetric Bernoulli random variables 
£™ :Q -> {-1, 1} with P(£ l = 1) = = -1). By Donsker's theorem, the sequence of 
random walks = ^= J2i=i £i"' > converges weakly to a standard Brownian motion 
B = (B t )te[o,i] [5], Theorem 16.1. The idea of Sottinen [21] is to combine these random 
walks with a pointwise approximation of the kernel in representation (15). Define the 
pointwisc approximation of z(t,s) as 

z^(t,s):=n [ S z(^-,u)du. 

Js-X/n V " / 

The sequence of binary random walks 

r t V nt \ r i/n 

B?> n := / zM(t,s)dB^=J2n 

Jo l=1 J(i-l)/ n 



in 



then converges weakly to a fractional Brownian motion te[o,i] m the Skorokhod 

space L>([0,1],R) [21], Theorem 1. 

A major advantage of the binary random walks is that we can avoid the difficult 
Wick product for random variables with continuous distributions. We approximate this 
operator on the binary random walks by discrete Wick products. 

For any neN, let (£",£2 > • ■ • >£n) be the n-tuple of independent and identically dis- 
tributed symmetric Bernoulli random variables for the binary random walk B^' n . The 
discrete Wick product is defined as 



ieAui 



ieA ieB I 0, otherwise, 
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where A,BC {1, . . . ,n}. We denote by 

the CT-field generated by the Bernoulli variables. Define 

He 

Clearly, the family of functions {S^ : A C {1, . . . , n}} is an orthonormal set in L 2 (n, T ni P). 
Since its cardinality is equal to the dimension of L (fl,J-' n ,P), it constitutes a basis. Thus, 
every X £ L (Q,J-' n , P) has a unique expansion, called the Walsh decomposition, 

AC{l,...,ra} 

where x\ € R. The Walsh decomposition can be regarded as a discrete version of the 
chaos expansion. By algebraic rules, one obtains, for X = J2ac{i n } x A^A an d Y = 

C{l,...,n} Vb^Bi 

Xo n Y= ]T ( E x aVb 

CC{l,..,n} \AUB=C 
AnB=<!) 

Furthermore, the L 2 -inner product can be computed in terms of the Walsh decomposition 

as 

E[XY}= J2 X AV% (17) 

iC{l...,Tl} 

There exists an analogous formula for the Wick product on the white noise space via chaos 
expansions that justifies the analogy between the discrete and ordinary Wick calculus 
(cf. Kuo [14]). For more information on the discrete Wick product, we refer to Holden et 
al. [10]. More generally, the introduction of a discrete Wiener chaos depends on the class 
of discrete random variables (£" , Q, . . . , £") . We refer to Gzyl [9] for a survey of other 
discrete Wiener chaos approaches. 
The representation 

Bf = £ bl& with bl % := ^ z(^-,s) ds 

~[ J(i-l)/n \ n J 

is the Walsh decomposition for the binary random walk approximating B H in 
L 2 (il, T n , P). Note that = 6™ nt j . i . Thus, we can consider B t ' n = B^l, as a pro- 
cess in discrete time. We can now state our first convergence result. 

Theorem 1. Suppose that: 
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1. lirrin^oo a n ,k = a-k exists for all k <G N; 

2. there exists a C G R+ such that |a„.fc| < C k for all n,k € N. 

The sequence of processes X]fe=o ^ 2 ^-{B H ' n ) <>n ' k then converges weakly to the Wick power 
series f^Lo ^{B H f k in the Skorokhod space D([0,1],R). 

The proof is given in Sections 5 and 6. 

Consider now the following recursive system of Wick difference equations: 

U* ■» = + kU^- n o n (Bffi - B*? iy J, U?> n = 1, = 0, (18) 

for all I = 1, . . . , n and k € N. This is the discrete counterpart of the recursive system of 
SDEs in (11). We observe that U°' n = 1 = (B H > n f^ and U hn = (B^")*" 1 , but 

rr2,n _ offri off,™ / r>H,n r,H,n / n ff,n\©„2 
2 1/n 11 2/n ' 2/n 11 2/n V 2/n ' 

Thus, in contrast to the continuous case in (11), the discrete Wick powers are not the 
solutions for (18) if k > 2. 

However, we can prove a variant of Theorem 1, based on the system of recursive Wick 
difference equations, whose proof will also be given in Sections 5 and 6. 

Theorem 2. Under the assumptions of Theorem 1, define U t ' n := as the piecewise 

constant interpolation of (18). 

The sequence of processes X)fc=o ^-jj L U k,n then converges weakly to the Wick power 
series J2h=o T\( BH ) <>k m the Skorokhod space D([0, 1],R). 

Example 1 (Wick powers of a fractional Brownian motion). For a n ,k = ^l{fc=z} ; 

U l ' n A {B H f. 

Example 2 ( Geometric fractional Brownian motion). For a n ,k = a k = 1, we have 

L«*J , 

exp »(B t ' ) := 2^ Y\ y 1 ' P *• 

k=0 

S n := T i U k ' n A exp°(B ff ). 

fe=0 

Observe that by summing up the recursive system of Wick difference equations (18), we 
obtain 

sr = sr_i + sr_ x <>,, (<„- - s? = i, (19) 
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for I = 1, . . . ,n, where SJ l = Sy n - Hence, the piecewise constant interpolation of (19) 
converges weakly to the solution of the fractional Doleans-Dade equation (13). 

The reasoning of the previous example can be generalized as follows. 

Theorem 3 (Linear SDE with drift). Suppose n,s 6t. a > 0. Then := SV, , 
where S n is the solution of the Wick difference equation 

S? = (l + SP-i + oSl-x <>n (<;" - B?^ 1)/n ), S^ = s , 1 = 1,..., n, (20) 

converges weakly to the solution of the linear SDE with drift 

dS t = / iStdt + o-S t d°Bf, S =s , (21) 
in the Skorokhod space Z?([0,1],R). 

Proof. First, observe that for o~ n — > a > and a nt k = cin^o'!^, we obtain, by Theorem 2, 
that 

n oo 
fc=0 k=0 

With the choice a n ,k ■= 1 and 

a 

<7„ := ; > a as n — > oo, 

1 + n/n 

we observe by (18) that V™ := V l r J n satisfies 

vr = VT-i + (tTJ^) ^ ° n {B "'" ~ B f-^ v ° = x - 1 = x > • ■ • - n - 

Consider now the piecewise constant function (W")t e [o,i] determined by W™ := Wy nt ^ 
and 

Wf=(l + ^Wl, ^o" = *o, Z = l,...,n. 
By this well-known Euler scheme, 

(W?) t6[ o,i] — >3 (exp(^)) te[0 ,i] (22) 
in the sup-norm on [0, 1]. The product 



[l + ^VT-lWT-y. 
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VfWS = s , 1 = 1, 



satisfies the Wick difference equation (20) for S*™ = V^W™. The multiplication by the 
deterministic function sq exp(/ii) is continuous on the Skorokhod space. Thus, with (22) 
and Billingsley [5], Theorem 4.1, we obtain 

0S?)te[o,i] = (V t n W t n ) te[ 0A] soiexp^exp^aB? )) tg[0il] 

in the Skorokhod space D([0,1],R). As s exp(/xt)exp°(CTBf ) solves the SDE (21) 
(cf. Mishura [15], Theorem 3.3.2), the proof is complete. □ 

Remark 1. Theorem 3 holds with additional approximations (cr n ,fj, n ) (a, //). 

Remark 2. Theorem 3 was conjectured by Bender and Elliott [3] in their study of the 
discrete Wick- fractional Black-Scholes market. They deduced an arbitrage in this model 
for sufficiently large n. Although the arbitrage or no-arbitrage property is not preserved 
by weak convergence, this model showed that it is even possible to obtain arbitrage in 
this simple discrete Wick fractional market model. In a recent work [23] , Valkeila shows 
that an alternative approximation to the exponential of a fractional Brownian motion by 
a superposition of some independent renewal reward processes leads to an arbitrage-free 
and complete model. We refer to Gaigalas and Kaj [8] for a general limit discussion for 
these superposition processes. 

Theorem 4 (Linear system of SDEs). The piecewise constant interpolation 

( vn v n \T ■ I yn V n S I T 

l A t i Y t ) ■— l A [ntJ > Y [nt\ ) 

for the solution of the linear system of Wick difference equations 
X? = X'f l _ 1 + {A X X?_ X + A 2 Y{L x )o n {Bff: - Bf'",\ X^ = x , 1 = 1,..., n, 



J l/n -"(i-l)/n/' 
J l/n _ a {l-l)/nh 



= Y?_ x + {B X X?_ X + B 2 YP_ X ) o n (B?£ - B^J, F " = y , l = l,...,n, 



converges weakly to the solution (X, Y) of the corresponding linear system of SDEs (1) 



in the Skorokhod space D([0, 1],R) 2 . 



Proof. Analogously to (14), we obtain, by the recursive system of Wick difference equa- 
tions for U k ' n in (18), the coefficients for the solution of the systems of difference equations 

OO OO j 

fe=0 k=0 
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recursively by 

a =x , b = y , a,k = Axdk-i + A 2 bk-i, h = B^au-i + B 2 bk-i- 
We define the upper bound 

AfAB:=2max{|Ai|,|A a |,|Bi|,|B3|}. 

Suppose that r\,r 2 € R arc arbitrary. By the linear system and (18), the sequence of 
processes 



k=0 

fulfils the conditions in Theorem 2 with 



\ria k + r 2 b k \ < max{|a;o|, M}(N + \r 2 \)M k AB . 



Thus, we obtain the weak convergence 



X n + r 2 Y n -A j[] ( riak t r ) = riX + r ^ Y - 



k=0 



The Cramer- Wold device (Billingsley [5], Theorem 7.7) can now be used to complete the 
proof. □ 

Remark 3. Theorem 4 can be extended to higher-dimensional linear cases. It also holds 
for an additional approximation of the coefficients A n ^ — > A4 and B n ^ — > Bi for n — > 00. 

Example 3 (Wick-sine and Wick-cosine). The piecewise constant interpolation of 

X? = Xf-i + Yf-i <>n (<r - B?»), X- = 0, l = l,...,n, 



J l/n "(/-l)/ 

converges weakly to the solution of the linear system 



Yr = Yll x Xf-i o n (Bffc B?*), y " = 1, / = 1, . . . , n, 



dX t = Y t d <> B?, X Q = Q, 

dY t = -X t d <> B?, Yb = l, 

the process (sm (B^ 1 ), cos°(i?f f )) T . By Theorem 1, it can also be approximated by the 
discrete Wick version functional (sin " (flf '"), cos " (B? '™)) T 



402 



C. Bender and P. Parczewski 



4. Walsh decompositions and Z 2 -estimates 

In this section, we give the Walsh decompositions for the approximating sequences and 
obtain some L 2 -estimates. A key to the approximation results will be the convergence of 
the L 2 -norms of the discrete Wick powers of Bf' n to the corresponding L 2 -norms of the 
Wick powers of Bf. 

Recall the Walsh decomposition B^' n = Xa=i ^t.i^i ■ Define 

h n — TT A™ «™ — TT t n rl n — h n h n 

U^11"U' ■— 11 Si' a l,i ■— °l/n,i ~~ °(l-l)/n.i 

for I = 1, . . . , n. Note that df^ = ft™/ v df^ = 0, for i > I and that the increment has the 
representation 

I 

r>H,n tjH,ii _ ST^ m tn 
i=l 

Recall the recursive system of Wick difference equations, 

E/f." = Ul\ + kU^ n o n - Bf» l)/n ), U?< n = 1, E#" = 0, (23) 

for 1=1, ... ,n and k £ N. 



Proposition 1. For all n, k £ N and I = 0, . . .,n, we have the Walsh decompositions 

J\ U l' n = E ( E II^H (24) 



CC{1, ...,1} x m:C-^{l,...,Z}peC 
\C\— k injective 



l. (B fp^ = £ b ^ cEh (25 ) 



CC{1,...,I) 
\C\=k 



liK:r nk -h u *' n= E ( E n (26) 

cc(i,...i) v m:C-s-{i, ...,;} pec 

|C|— ft not injective 



Proof. We use the conventions that an empty sum is zero, an empty product is one 
and that there exists exactly one map from the empty set to an arbitrary set. For these 
reasons, the formulas hold for k = or I = 0. We prove (24) by induction as follows. For 
alU = 0, . . . , n and all k £ N, it is obvious that J7 ; '™ = 1 and U ' n = 0, as in formula (24). 
Suppose the formula is proved for all positive integers less than or equal to a certain k 
and all I = 0, . . . , n. Furthermore, for k + 1, suppose the formula is proved for all positive 
integers less than or equal to a certain /. For k+ 1 and I + 1, we compute, by the difference 
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equation (23) and the induction hypothesis, 

u^n_ U k + i, n = {k+1)k J £ £ n^),, s 'c)o«E d r +M er 

Vc{l,..J} m:C-^{l,...,/}peC ' i=l 

|C|=fe injcctivc 

= (fc + l)! £ XI II CwX+LiSSuW ( 2? ) 
cc{i, ...,;} • .{ : i\ /r < • 

ig{l,...,i+l} injective 
\C\=k,i(£C 

=(*+i)i e n^),^- 

C'C{l,...,i+l} m':C'-»{l,...,J+l} pec 
|C"|— injcctivc, 3q:m(g)— Z + l 

Note that d m ^ p = for all p — 1 > to. Thus, by the induction hypothesis, 

c/, fe+i '" = (fc+i)! E IPmw^s- ( 28 ) 

cc{i,...,i+i} ».:( • -i i / • 1 1 pec 

|C| = fc+l injcctivc, Vq:m(q)<l+1 

Thanks to equations (27) and (28), wc obtain (24). In particular, U,' n — for all k > I. 
We now compute the fcth Wick power of (B l as follows: 

/ n \ °" k i / k k \ 



£w? = E n *?/».<, n 3 

pairwisc distinct 

e wfn^nff)= e 

{!,...,/} \i6C iec / CC{1,.. 



CC{1, ...,/} VGC ieC / CC{l,...,i} 

|C|=fc |C|=fc 

In particular, {Bf.^f^ = for all k > I. This yields (25). 
The telescoping sum yields 

i i i 

E/ ^™>(p),P = E/ ^»™(p),P = E/ (^m(p)/ra,p — b(m(p)-l)/n,p) ~ ^l/n,p 
m (p) — 1 m(p)—p m(p)—p 

and thus we get 

e n <(p),p = n ( e c ( p),p ) = n ^ = ^.o- (29) 

m:C->{l,....l}pGC pSC \m(p) = l / peC 

Equation (26) is, thus, implied by (24) and (25). □ 
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In the following propositions, we obtain some elementary estimates for the L 2 -norm of 
discrete Wick powers of Bf' n . 

Proposition 2. For all t G [0, 1] and i £ {1, . . . , [nt\ }, 

bl l <2c H n- {1 - H) . 

Proof. We estimate 

V 2 / Jfi-l)/n Js 



<n 1/2 c H ( H- 



< n 1/2 c H 



(i-l)/ 
pifn 



H ( l nt \ 



1 



3/2 -H 



(i-l)/n 

3/2-H 



H-l/2 



1 / [nt\ 



H-l/2\ n 
l\ 3/2 - H \/[nt\\ 2(H - 1/2) 



H-l/2 



ds 



1 



Since t < 1, jj^jj < 2 and \x\ 3 ^ 2 - H - \y\ 3/2 ~ H <\x- y\ 3/2 ~ H , the assertion follows. □ 



Remark 4- Observe that 



E[flf ,n Bf = E 



.11 ,12 = 1 



LntJ 



(30) 



By Nieminen [19], we thus get, for any s,i£ [0, 1], the following convergence: 



IMI r i/r 
nB H.n B H, n] = / 

i=l J ( i ~ 



/(i-l)/r 



[nt\ 



du 



i/n 



(i-l)A 



[nsj 



— > [ z(t,u)z{s,u)du = E{Bt I B?}. 
Jo 

Moreover, we have, by the Cauchy-Schwarz inequality, the upper bound 



(31) 



l nt l / r i/r, 

E[(B^-B^) 2 ] = J2[V^ 



l nt l r i/n 



(i-l)/n\ V " 

LntJ 



u — z 



du — z 



[nsj 



[nsj 



dt/ 



d(/ 



(32) 



n n 



2H 
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Proposition 3. For all t> s in [0, 1] and all N EN such that [ns\ > N , we have 



< E[(5f > n f] N +E[(Sf ' n f] N - 2E[B?> n B?< n } N 



(33) 



<2c 2 H N 2 t 2H( - N -^n- (2 - 2H l 



In particular, 



H,n\o n N _ / r>H,n\o n N\2 



lim J&[((B?> n y»«-{B 



} = n((B?y N -(B?y N f}. 



Proof. As N < [ris\ , we get, making use of Proposition 1. (17) and (30) in Remark 4, 



A! 



E[(Bf'") " JV (Bf' n )°" JV ] 



AH 



-E 



CC{l,...,|_nt|} ' ^ CC{l,...,|nsJ} 

|C|=JV |C|=7V 



s,C 



CC{l,...,|ntJ} 
|C|=JV 



[_nt\ N Y nt \ N 

= e k< • e n^- 6 ^) 

ii,...,ijv = l j = l ii,...,ijv = l J=l 

\nt\ N 

= ElB?' n B^T- e n^w 

ii,...,ijv = l J=l 
3fc,Z : ik=il 



Thus, we have 



(34) 



E[(B"' n ) 2 ] N +~E[(B"' n )' 2 ] N - 2E[B^ n B^ n ] N 

-±E[((B^r^-(B^r- N ) 2 } 



N 



N 



\nt\ / N 

E II '< " • 1R 2 1R- K, ^ 

i u ...,i N =l \j=l j = l j = l 



(35) 
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[nt] / N N \ 2 

= e nw-nw*) 

it,...,i N = l \j=l 3=1 / 

: ik—ii 

Hence, the left-hand side of the inequality in (33) follows. By Proposition 2, (30) and 
(32) in Remark 4, as well as < t, we obtain 

[nt] 



[nt] / N N \ ■ 

e iiw-iiwa) 

..,i N = l \j=l j=l I 



3fc,Z : ik=ii 



N 



[nt] 

* E (Efe) 

ii,...,ijv=l V/=l / 
3fc,i : ifc=ii 



< 



N 2 J (max(6^) 2 



[ntj /jv-1 \ : 

e n&*> 

.,ijv-i=l \ j=l / 



(36) 



< 2(^JV 3 E[(Sf '")2]W-l n -(2-2H) < 2c 2 ?iv 2 t 2H(iV-l) n -(2-2H) ^ g 

for n — > oo. The representation of Wick powers of i?^ by Hcrmite polynomials, as in (9), 
their orthonormality (cf. Kuo [14], page 355) and the polarization identity collectively 
yield E[(Sf ) oAr (Bf ) oJV ] = iV!E[(Bf )(Sf )] w (cf. also [20], Lemma 1.1.1). Thus, we have, 
by (35), 



mB H,ny n N_ [B H 



2 }-n((B?y N -(B»y N y 



N\(E[(B?' n f] N - E[(5f ) 2 ] w + E[(Bf •") 2 ] JV - E[(Bf ) 2 ] w 



- 2E[B?' n B? - n ] A + 2E[Bf Bf ]*) 
L«*J / jv jv \ 2 

e n(^)-iK^) ■ 

ii,...,ijv=l \j=l J=l / 

3fc,i : ifc=i( 

Applying the convergences (36) and (31) yields 

E[((Bf - (Sf ' n )*" JV ) 2 ] - E[((Sf ) oAr - (Bf ) o7V ) 2 ] 
Remark 5. In particular, we obtain, by (34), (32) and |-^^-| < t, 

E (&"c) 2 = (^ I ) 2 E[((^"r" Ar ) 2 ] 



0. 



CC{l,...,|nt|} 
|C|=JV 



<^E [( ^)r<± 



2ffN 



2HN 



□ 



(37) 
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The next proposition estimates the difference between the approximating sequences in 
Theorems 1 and 2. 



Proposition 4. Under the assumptions of Theorem 1, there exists a constant K > 
such that for all t € [0, 1], n > 1 and k € N, 



\ " a n,k i p _H\n^o„fc \ ~* a n,k f T k,- 
k=0 k=0 

for the approximating processes in Theorems 1 and 2. 



(38) 



Proof. Recall that d^ = b« /n . - 6f r _ 1)/n> . = V^I^yM^ s ) ~ <^s)) ds. By (32) 
in Remark 4, we obtain 



«j 2 <E«^ 2 ^ 



r r — 1 



2fl 



-2H 



Thus, we have d r ^. i <n H for all i,n,r> 1. Hence, we obtain, as the sum in (29) tele- 
scopes, for \C\ > 2, 



e 

m:C^{X,...,\nt\}leC 
not injective 



E II d ™(0,; 

m:C->{l,...,|ritJ} Z6C 
3u,vGC : m(u)— m(f;) 



E E 



n < 



,(/),; 



iiGCm:C\{u}-*-{l,...,|ntJ} x !eC\{«} 



E Cwi 39 ) 

i>ec\{u} 



<max^(|C|-l) ^ 



E II 

C'CC m:C"->{l,...,[niJ}ZGC" 
|C'| = |C|-1 



<n- H (\C\-l) ]T b lc>- 



C'CC 
|C'| = |C|-1 



By (26), (39), (37) and since ([nt\ — (k — 1)) < n, we obtain, for k> 1, 



E 



/V r>H,n\o n k Tr k,n\ 
) ~ U t ) 



< 



CC 



E (n- H (k-l) £ 

{!,..., L«*J } C'CC ' 



\C\=k 



|C'|=|C|-1 
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<n-™(k-l? Y, ( fc -!) E Kef 

CC{1 [nt]} C'CC 

\C\=k |C'|=|C|-1 



< n™(k - lf([nt\ -(k- 1)) ( & "c<) 2 < Jj-J 

C'C{l,...,LntJ} 
|C"|=fc-l 



t 2H(k-l) n l-2H 



Since a nj fc < C k in Theorems 1 and 2, and ((Bf ' n ) ™ fc — U^' n ) are zero for k = 0, 1 and 
orthogonal for different k by Proposition 1 and (17), we get 



E 



E a n,k , u H,n^ n k a n,k f%k, 



k=0 
fc=2 



fc=0 



a n,k /Y D-H,n\o„fe [ft; 



2n 



< 



E^ 

\fc=2 



2/,- 



(fc-1)! 



2H(k-l) „l-2ff 



As the series on the right-hand side converges uniformly in t £ [0,1], the assertion fol- 
lows. □ 



5. Convergence of the finite-dimensional distributions 

We first prove that Theorems 1 and 2 hold with weak convergence replaced by convergence 
of the finite-dimensional distributions. To this end, we first approximate the Wick powers 
of Bf by induction. We then combine these convergence results to approximate the Wick 
analytic functionals F°(Bf ) = Ysk=o TfeK-^t ) • Finally, we conclude that convergence 
in finite dimensions holds in Theorem 2. 

We observed in Section 2 that (Bf ) oJV = /ij^ 2Jf (Bf ) and that the Hermite recursion 
formula 

{B Hy(N+l) = {B H ){B HyN _ \ t \™ N {B Hy(N -1) (4Q) 

holds. For the discrete Wick powers of the discrete variables, we now obtain a discrete 
variant of (40). 

Proposition 5 (Discrete Hermite recursion). For all N > 1 and t e [0, 1] , 

(Bf = B?' n (B?> n y« N - NE[(B?> n ) 2 ](B?> n y^ N -V +R(B?' n ,N), (41) 

with remainder 

R(B?> n ,N) = N\ ]T &£c3&E( 6 m) 2 ( 42 ) 

CC{l,...|ntJ} iSC 
\C\=N-1 
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and 

E[(R(B? AO) 2 ] < lQc%N\N 3 n- ( - 4 - 4H \ (43) 

In particular, we will use the fact that the discrete Hermite recursion (41) converges 
weakly to Hermite recursion (40) for n — > oo. 

Proof. By Proposition 1, we get 

{B H, nyn{N+1) = B H,n ^ {B H,ny nN = fj^ ^ ^ ^ f £ Nlb^) 

\i=l / MC{l,...,|ntJ} ' 

\A\=N 

(44) 

= Bf >»(Bf -EE mh ii h umt- 

A<Z{l,... 1 \nt\}i&A 
\A\=N 

For the second term in the last line in equation (44), by (30) in Remark 4 and Proposition 
1, we obtain 

E J2 mb li b um? 

A<Z{l,...,\nt\} ieA 
\A\=N 

= m E E 6 U« S W & "^) 2 ^ ! E Kc^cT,^) 2 

AC{1 lnt\}ieA CC{l,...,|ntJ} i£C 

\A\=N |C|=JV-1 

/|nf| 



/[nt] \ 

N(N-1)\ E b lc^c E( 6 m) 2 "E( 6 m) 2 
cc{i,...,|ntj} \i=i iec / 



2 



|C?|=JV— 1 

= AT(Bf '")*-^-i)E[(Bf '") 2 ] — AH E b tc^c E( 6 " 

CC{l....,|ratJ} iGC 
|C|=iV-l 

which yields (41) and (42). Thus, thanks to Proposition 2, Remark 5 and t < 1, we obtain 



2 

2 ' 



E[(i?(i?f AO) 2 ] = (TV!) 2 E (Kc) 2 E( 6 m) 

CC{l,...,|ntJ} MSG 
|C|=iV-l 

- (7V!)2 (jv^T)T* 2H(JV " 1)((JV - i) 4 ^^ 2 ^) 2 

< 164A^!Ar 3 n-( 4 - 4ff ). □ 
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Theorem 5. For all JVeN, 

(l,B H ' n , (B H ' n ) <>nN ) A (1, B H , . . . , (B H )° N ). 



(45) 



Proof. The proof proceeds by induction on N. By Sottinen's approximation, (1, B H ' n ) — > 
(1,B H ). Suppose that equation (45) is proved for some JV > 1. Assume that k <G N and 
r\ £ M for j = 0, . . . , AT + 1 , i = l,...,k and ii, *2, . . . , tfc S [0, 1] are chosen arbitrarily. By 
the pointwise convergence ~E[(B^' n ) 2 ] — > |i| 2ff and the generalized continuous mapping 
theorem (Billingslcy [5], Theorem 5.5), the induction hypothesis implies that 



N / k 
i=0 \j = l 



o,,/ 



^rf+ 1 (<'"(<-)^ w -iVE[(<») 2 ](<' n ) " (JV ~ 1) ) 



N / k 



1=0 \j=l / j=l 

Since H > §, (43) yields R(B^ n ,N) — > in L 2 {n,P). Thus, by Slutsky's theorem [5], 
Theorem 4.1, and the Hermite recursions (40) and (41), we obtain 



JV+l / k 



JV+l / k 



oZ 



z=o \j=i / z=o \j=i / 

By the Cramer-Wold device (Billingslcy [5], Theorem 7.7), we have 

(l,B H ' n , . . . , (B».n)*-JV+i) A (l, b h , . . . , (B H ) oN+1 ) 
and the induction is complete. 



□ 



Proposition 6. In the context of Theorem 1, convergence holds in finite- dimensional 
distributions. 

Proof. By Billingslcy [5], Theorem 4.2, it suffices to show that the following three con- 
ditions hold: 



fe=0 ' A;=0 



k \ v-t / ■ fc! v asn^oo; 



Vte [0,1], lim limsupE 



fe=0 fc=0 



Al 



(46) 
0; (47) 
(48) 
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Condition (46) follows directly from Theorem 5 and the generalized continuous mapping 
theorem ([5], Theorem 5.5). For the second condition, we compute 



E 



E a n,k , „H,nsok a n,k , D _f/,raxc 



\k=Q 



E 



. \fc=m+l 



fe=0 

Q-n.fe / T)H,n\Q n k 



E 



E 

fc— m+l 



~fcT 



E[((Bf.»)*-*)2]< £ 



fc!f 



applying the estimate of Remark 5. Here, we used the fact that discrete Wick powers of 

i+i ~fcT = 



different orders are orthogonal. Thus, we even obtain linim^oo limsup^^^ Y]^_ m , x - 



and, for all t G [0, 1], a stronger result than condition (47), 



lim lim sup E 



fci 



L \fc=0 fc=o 
By the orthogonality of the different Wick powers, we have 
\ 2 

a k i j^H \ok 

\k—m-\-l 

for m— ?> oo, which implies that condition (48) even holds in L 2 (il,P) 



= 0. 



E 



E Sr^ 



oo / \ 2 

afc 



E 

A;— m+l 



fc! 



E[((i?f )♦*)*]< £ 

/c— m+l 



2A- 



c 
~k\ 



,2Hk 



□ 



In view of Proposition 4 and Slutsky's theorem, the previous proposition also implies 
the following. 

Proposition 7. In the context of Theorem 2, convergence holds in finite- dimensional 
distributions. 



6. Tightness 

We now show the tightness of the sequences in Theorems 1 and 2 by the following 
criterion, which is a variant of Theorem 15.6 in Billingsley [5]. 

Theorem 6. Suppose that, for the random elements Y n in the Skorokhod space 
D([0, 1],K) andJ2Z ^(B H r k m C([0,1],R), 

oo 

yn f d > a k ^Hyk 
k=0 ^' 
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and for s <t in [0, 1], 



m[(Y?-Y?y]<L 



[nt\ \ns\ 



211 



where L > is a constant. Then Y n converges weakly to X^fcLo it (B ) in D([0, 1], 
Proof. Let s < t < u in [0,1]. By the Cauchy-Schwarz inequality, 

E[|y t "-y»irc-y»|] 
< (E[|y t " - r;f ]) 1 / 2 (E[|y„- - y/T]) 1/2 



<L 



\nt\ \ns\ 



n 



\nu\ \nt\ 



i) 



n 



H 



[nu\ [ns\ 



i) 



n n 

If u — s > i we have, since [nu\ < nu and — \ns\ < —ns + 1, 

2ff 



2 if 



\nu\ \ns\ 



n 



n 



<(2(u-*)) 
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and thus 



E[|y t n - y; i ||y„ n - *?0 < - a) 



\2H 



(49) 



If u — s < — , then we have either [nsj = [nt\ or [nij = \ nu\ and so the left-hand side in 
(49) is zero. Thus, the inequality (49) holds for all s < t < u. By the convergence of the 
finite-dimensional distributions and [5], Theorem 15.6, we get the weak convergence of 
the processes. □ 

For the application of this criterion to the discrete Wick powers, we need two lemmas. 

Lemma 1. Let [X, (-, ■)) be a real inner product space and \\x\\ :— (x, x) the corresponding 
norm on X . Then, for all x, y G X and N > 1, 

\\x\r + \\y\r 2((x,y)) N < 2 N+ \\\x\\ + \\y\\f(^)\\ x - yf. 
Proof. It holds that 



>,[\(\\x\\ 2 + \\y\\ 2 -\\x-y\\ 2 ) 



N 



2 N-1 



N-l , s 

(IMI 2 + \\y\\ 2 ) N + E U dl-ll 2 + \\y\\ 2 ) k (-i) N - k \\x y\\ 2iN - k) 

k=0 ^ ' 
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HI*-2/II 2 ^ztX;( ^(INP + lbll 2 )^-!)^- 1 ^-^! 2 ^- 13 - 

fc=0 ^ ' 

Hence, we get 

\\xf N + \\vf N -^y)) N 
= M 2N + \\y\\ 2N -^T(M 2 + \\y\\ 2 ) N (so) 

fe=0 ^ ' 



Since (^) N 1 J2k=o (k) = ^> * ne nrs * nne on * ne right-hand side of (50) can be treated 
as follows: 

INI 2JV + NI 2JV -^T(lkll 2 + IMlT 

(51) 



JV-l 



2fc|, ,||2(iV-fe) 



2 N-1 \ k 

k=l v 

As ( , ) = \ N _u), we now collect the corresponding summands in sum (51) for f . We 
obtain, by the mean value theorem with 

M x , y := max (A||s|| + (1 - A)||y||) = max^U, ||y||} 



Ae[o,i] 



and since k(N — k) < 



4 ' 



lkll 2Ar + l|y|| 2Ar -NII 2fe bll 2(7V - fc) -lkll 2(Ar - fc) !l2/ll 2fc 
= (bll 2fe -!|y|| 2 ' ; )(bll 2(Ar - fe) -||2/ll 2(JV " fc) ) 

< 2kM* k - 1 \\x y\\2(N- k)M^ y N - k ^\\x y\\ < N 2 M^ y N ~^ \\x y\\ 2 . 
Analogously, we obtain, for k = ^ , 



M2N + M * N \\x\f\\y\\ N = k\\x\\ N \\y\\ N f < \m^-Vn*\\x 



N + h\\ 2N 

2 311 2 y ' ~ 2" 

Plugging these estimates into (51) and recalling that (\) N_1 J2k=i = 0- ~ pr=r)j 

we obtain 

Nr + \\y\\ 2N \\\*\\ a + h\\ 2 ) N < (i - ^)n 2 mW-V\\x- y || 2 . (52) 
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For the term in the second line on the right-hand side of (50), we observe that, by the 
triangle inequality, 

(\\x\\ 2 + \\y\\ 2 ) k (-^ N - k - 1 \\^-y\\ 2{N - k - 1) 
< (11-11 + IMI) 2fc (NI + \\y\\) 2{N - k ^ = (H + IN) 2 ^. 



Applying 

M x , v < \\x\\ +II2/H, 

and (52) to (50), we have 

\\x\\ 2N + \\y\\ 2N -2{{x,y)) N 



N-l N-l 



E 

k=0 



JV 



= 2- 



1 



2 N-1> 



< 



1 - 



2 N-1 



JV + 2 — 



2 N-1 



(\\x\\ + \\y\\)^ N ^\\x-y\\ 



By a short calculation and induction, we obtain 

2 ~ ) " 1 {^3}2 Ar + 1 { jv=3}^ < 2 W+1 . 



1 - 



2 n-i 



N 



□ 



Lemma 2. i*br all t> s in [0, 1], we /iaue 



JV! 



E[((Bf - (flf • n )°" N ) 2 } < 8 N 



\nt\ \ns\ 



n n 



2H 



Proof. For JV = 1, the inequality is fulfilled by (32) in Remark 4. For JV > 1, we consider 
the cases iV > [wsJ and [nsj > JV separately. For |_ri£j > JV > [nsj , we have (^B^' n ) <>nN = 
0. Hence, Proposition 3 and Remark 5 imply that 



1 

JV! 



Since JV > 2, 2H > 1 and < 1, wc obtain 



2HN 



[nt\ 



2HN 



< 



[nt\ 



(([nt\-[ns\) + [ns\Y 



< ((L"tJ-Ln*J) + ^) < " ((Ln*J-L"*J)(^ + l)) 



(JV+1) 2 



[nJj L ns J 



n n 
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Since (N + l) 2 < 3^ for N > 2 and 2H < 2, we obtain 



N 



[nt\ \ns\ 



2H 



for all \_nt\ > N > [ns\ . Recall that, by Proposition 3, for \nt\ > [ns\ > N, 

±E[((Bf' n r" N -(B^ n r N ) 2 } 

< E[(Bf ' n ) 2 ] N + E[(B^ n ) 2 ] N - 2E[(Sf ' n )(Bf > n )] N . 
For any n £ N, we can rewrite 



415 



E[(Bf'")(Bf>™)]=£&^ 



i=i 



as an ordinary inner product on R" of the vectors (bf 1 , . . . , 6"„) T and (6™ lf . . . , &™ ra ) 
Thus, the application of Lemma 1 with t,s£ [0, 1] gives 



1 

m 



E[((#f ■ n )°" Ar - (flf '")°" JV ) 2 ] 



< 2 JV+1 (E[( J Bf' n ) 2 ]V 2 +E[(Bf'™) 2 ]V 2 ) 2 ( JV - 1 )E[((Bf' n ) - (Sf™)) 2 ] 

< 2A r +i2 2 ( JV - 1 ) 



1) 


[_ntj [_n-s\ 


2H 

< 




\nt\ \ns\ 


2H 




n n 






n n 





If N > \ nt\ , then the left-hand side of the assertion vanishes. 

Remark 6. The proofs for a fractional Brownian motion on some interval [0,T] C 
follow by a straightforward modification. As E[(fif '") 2 ] < T 2H for i G [0,T] and 



□ 



2HN 



< rp2H(N-l) 



[nt\ 



2H 



we obtain the previous lemma for t > s in [0,T] as 



AT 



|_nij [ ns J 



2// 



We are now able to prove the weak convergence to the Wick analytic functionals of a 
fractional Brownian motion. 

Proof of Theorem 1. We apply Theorem 6. The convergence of finite-dimensional 
distributions was shown in Proposition 6. Let s < t in [0, 1]. Recall a n ^ < C ■ Then, by 
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the orthogonality of ((B^ I ' n ) <>nk — (^B^' n ) <>nk ) for different k and Lemma 2, we have 

\ 2 n 



E 



a n,k , D H,nxo n k a n,k , u H,n\O n k 

\k=0 k=0 j 



k=0 



= E( TT ) E[(«'T" fe - (flf "T" fe ) 2 ] < E ? 

fc=0 ^ ' ' 

Since < ^^L ^~W~ = exp(8C 2 ) =: L < oo, we have 



\nt\ \ns\ 



n n 



211 



E 



a n,k i D _ff.riNO„fe a n,k , u H,n\o n k 



fc=0 



r 


< L 


\nt\ \ns\ 






n n 



2tl 



(53) 

□ 



The alternative approximation, stated in Theorem 2, follows similarly, as we shall now 



sec. 



Proof of Theorem 2. Let s < t in [0, 1]. Recall that d 7 ^ i > only if i < m. Thus, by 
Proposition 1, we can write 



£W n =£«M £ E IP 



|ntj 



/v 



k=0 C<Z{l,...,[nt]} m:C-t{l,...,[ns]} leC 

\C\=k injective 



Observe that, by the telescoping sum in (29), we have 



e n< w> *< 



m:C-t{l,...,lnt\} l€C 
injective 
3u:m(u)> [nsj 



e n d ™(i),i 

m:C->{l y n t\}iec 

3u : m(u)> [ ns J 

= e n d m(o,i _ e n d ™(i),i 

m:C^{l lnt\}leC m:C->{l,..., [n«J } ieC 

Thus, due to the orthogonality of U t ' n — U k,n for different values of k, Proposition 1 and 
estimate (53), we obtain 

/ n n \ 

Ett^-Ett^" 



E 



\k=0 



k=Q 
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\nt\ 



E<* E E net,,,, 



fc=0 CC{l,...,|ntJ} v m:C^{l,...,LniJ} ItC 

\C\—k injective 

3it:m(ii)> [nsj 



[nt] 

<E<* E ( fo "c-^c) 2 <^ 

fc=0 Cc{l,...,|nt|} 
|C|=fc 



2/f 



and the result follows from Proposition 7 and Theorem 6. 



□ 
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